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The stabihty analysis of an anisotropic inflationary universe of the four dimensional Neveu- 
Schwarz-Neveu-Schwarz string model with a nonvanishing cosmological constant is discussed in 
this paper. The accelerating expansion solution found earlier is shown to be stable against the 
' perturbations with respect to the dilaton and axion fields once the dilaton field falls close to the 

' local minimum of the symmetry-breaking potential. This indicates that the Bianchi I space tends 

04 . to evolve to an isotropic flat Friedmann-Robertson- Walker space. This expanding solution is also 

shown to be stable against the perturbation with respect to anisotropic spatial directions. 

^ ■ I. INTRODUCTION 

(N 

fSj , Observations of the cosmic microwave background radiation, galaxies and other astronomical objects reveals that 
J> our universe, on a very large scale, is remarkably uniform and is currently under accelerated expansion [Q-Q. 

Therefore, a physical acceptable cosmology should provide a self-contained mechanism to smear out the primordial 
' anisotropy and achieve an accelerating expansion at present time. Furthermore, it is firmly believed that the Einstein's 
, theory of general relativity breaks down at high enough energy which, at least, occurred during the early epoch of our 
' universe. It is known that the string theories are promising candidates for all particles and interactions, including the 
T-H _ gravity, in a unified formulation. Thus the astrophysical and cosmological implications of string theories have become 
■ a highly developing research themes (see |^ and references therein). 

\ An anisotropic cosmology was considered recently in the framework of four-dimensional Neveu-Schwarz-Neveu- 
"'Y' ' Schwarz (NS-NS) effective string theory in which the gravitational field is coupled with the dilaton and axion fields 
@- In the de Sitter geometry configuration, i.e. with a positive cosmological constant A > 0, an infiationary solution 
' can be found in an exact parametric form. At large time limit, this universe is isotropilized and it's expansion is 
accelerating which is consistent with our current astronomical observations. Stability analysis of any solution may 
provide more information of the theory under consideration pO|-p^. Therefore it is interesting, and is the purpose of 
the present paper, to investigate the stability conditions of this anisotropic inflationary string cosmology. Our analysis 
indicates that the inflationary solution found in Rcf. j^] remains stable when the potential U{4>) is close to the local 
minima of the potential. 

This paper is organized as follows. In Sec. II we briefly review the NS-NS effective string theory and the exact 
solution found in |^ . The stability analysis with respect to the perturbations of the dilaton and axion flelds will be 
performed in Sec. III. The perturbation with respect to the gravitational field will be studied in Sec. IV. In Sec. V, 
we will draw some conclusions. 



X 



II. ANISOTROPIC SOLUTION 



The four-dimensional NS-NS effective action, which is common to both the heterotic and the type II string theories 
[§, is given by 
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S = I d'x^g \ R - K{dct>f - ^e-'^ iJ^ -U\, (1) 



where (f) is the so-cahed dilaton and H^^x — d^^B^x^ is a totaUy anti-symmetric tensor. Moreover, the k denotes the 
generahzed dilaton couphng constant {k — 2 for typical superstring theories) and U = U{(j)) is a dilaton potential. 
The field equations of the action (pi) can be shown to be 



(2) 



\/^{e-^^H^"'^) =0, (3) 
W^<i^+^e-^*H^~^d^U = G. (4) 

In addition, the totally anti-symmetric tensor H obeys the Bianchi identity d[^H^xp] = 0. 

In four dimensions, any three-form can be mapped one-to-one onto its dual one-form. This mapping relates the 
three rank tensor to a pseudo-vector A via H = *A with A — Af^dx^ the dual one-form. Moreover, one can show 
that the tensor H can be solved to give 

H^""^ = e^-^ e'^"^" dph (5) 

following the field equation (||). Here the totally anti-symmetric tensor e denotes the Levi-Civita tensor and h — h{t) 
is known as the Kalb-Ramond axion. Note that we use the convention that e^^'^'^P = —5'^Q5\52 5'^^/\/—g in this paper. 
In addition, the above solution holds if and only if the first homology group for the space considered is trivial [p^|-pl| . 
Consequently, the Bianchi identity shown earlier becomes 

d^{^ge^^d^'h)=0. (6) 

We will consider the case that the potential of the dilaton field is a positive constant 

t/(</)) = A > 0, (7) 

which represents a de Sitter space. An exact solution for an anisotropic (Bianchi type I) cosmology has been found 
in Reference Note that the line element of a Bianchi type I space, an anisotropic generalization of the flat 
Friedmann-Robertson- Walker (FRW) geometry, can be written as 

ds^ = -dt^ + a\{t)dx^ + al{t)dy^ + al{t)dz^ , (8) 

with ai{t),i — 1,2,3 the expansion factors on each spatial directions. With the identities (||, 0), the field equations 
(^, ^, ^) will take the following form 

^ 1 1 

3^ + ^ ef + + Ae4^ 1^^-^^ = 0' (9) 

4=1 

^dtiVe,)-^A^O, z= 1,2,3, (10) 

h + 3eh + 4:^h = 0, (11) 

ldt{V^)--e'^h' = 0. (12) 

V K 

Here we have introduced the volume scale factor, V :— Y[^=i directional Hubble factors, Oi :— dij ai, i — 1,2, 3, and 
the mean Hubble factor, 9 := J^Ui ^i/S = V/3V for convenience. In addition, we will also introduce two basic physical 
observational quantities in cosmology: the mean anisotropy parameter, A := ~ 9)^/39^, and the deceleration 

parameter, q dt9^^ — 1. Note that ^ = for an isotropic expansion. Moreover, the sign of the deceleration 
parameter indicates how the Universe expands. Indeed, a positive sign corresponds to "standard" decelerating models 
whereas a negative sign indicates an accelerating expansion. 

It was shown that the general exact solution of this model can be parameterized in the following form j|] 
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a,(T) =a,osinh".^ ^cosh".^ ^, i = 1,2,3, (13) 
e2^(^) = ^g(tanrI+tanh--0, (14) 
/ s KV(lh) tanh^'^ i - 1 . . 

''M = ''"-^T^s;;;i?4TT- '''' 

Therefore, one has 

0{t) = y^cothr, (16) 

y(T) = Vbsinhr, (17) 

where T ^3A/2(t-to), := l/3±y/2/3AK^/Vo, OJ := \/8</>o/3A/Vb and := y'C' V8«:0o- Hereio, aio,ifi, C > 

0, (j>o > and ft-o are free parameters representing the constants of integration. In addition, one writes Vq = 

for convenience. One can also show that there is an additional constraint Ki — following the field equations. 
These integration constants also obey the consistency condition, 

3 

■.= Y,K!^AV,f~K^o, (18) 
1=1 

which follows from equation (^). Thereby the mean anisotropy and the deceleration parameter are given by 

Mr) = sechV, (19) 
g(T) =3scchV-l. (20) 

This exact solution indicates that the evolution of the Bianchi type I universe starts from a singular state, but with 
finite values of the mean anisotropy and deceleration parameters. In the large time limit the mean anisotropy tends 
to zero, A 0, and the universe approaches an isotropic inflationary de Sitter phase with a negative deceleration 
parameter, (j < 0, providing an accelerating expanding universe consistent with the present observations. Furthermore, 
in the large time limit the dilaton and axion fields become constants, limt^cxD ^(i) = = constant, limj^oo e^'^*-'^ = 
2(/5g = constant. Note that the dynamics and evolution of the universe is determined only by the presence of a 
cosmological constant (or a dilaton field potential). There is no direct coupling between the metric and the dilaton 
and axion fields other than the constraint on the integration constants. 

Note that the mixmaster model ||2^-^ discussed earlier deals with a model with perfect fluid. Earlier work focuses 
on the behavior of the solution near the singularity. The static field equation studied in Ref. |^ is similar to the 
mixmaster model with axion and dilaton behaving like the perfect fluid. The solution shown in this section is, 
however, the first time such exact solution is found for the NS-NS string effective theory that exhibits the desired 
property driving an initially anisotropic Bianchi type I space to an isotropic space. We will focus on the stability 
property of this solution in the following section. 

III. PERTURBATION AND STABILITY 

Perturbations of the fields of a gravitational system against the background evolutionary solution should be checked 
to ensure the stability of the exact or approximated background solution. In principle, the stability analysis should 
be performed against the perturbations of all possible fields in all possible manners subject to the field equations 
and boundary conditions of the system. In the following section, we will divide the perturbations into two disjoint 
classes: (a) the perturbations of the scale factors, or equivalently the metric field, and the axion field; and (b) the 
perturbations of the dilaton and axion fields. 

Note further that the axion field can be solved as a combination of the dilaton field and the scale factors according 
to Eq. dnl). The result is 



/i = Ce"'^'^y-\ (21) 

with a constant of integration C > 0. Here C — Q means that the axion field does not couple to any other field 
Therefore, the effect of the perturbation of the axion field can be replaced by the equation (pl[). In fact, Eq. (|2~ 
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indicates that if one perturbs the axion field, one has to perturb either the dilaton field or a; altogether unless the 
constant C vanishes. Note also that equation ( pT| ) indicates that h ~ which is much smaller than a in the large 
time limit where Oi ^ a for all directions. This is an indication that the effect of the axion field perturbation can be 
ignored in the metric perturbation we will discuss later. 

Once we replace the effect of the axion field in the field equations (^), (|lO|), and (|T^), one can further simplified 
the perturbations into two disjoint class: (a) the perturbations of the scale factors and (b) the perturbations of the 
dilaton fields. We will argue that the most complete stability conditions we are looking for can be obtained from class 
(a) and class (b) perturbations even the backreaction of the scalar field perturbation on the metric field perturbations 
is known to be important j27| . We will show that this backreaction does not bring in any further restriction on the 
stability conditions. 

The reason is rather straightforward. One can write the linearized perturbation equation as 

DlM,+Dl5c^^Q (22) 

for the system we are interested. Here the axion perturbations have been replaced according to Eq. Mor cover, 

perturbations are defined as 0^ = 0° + 5ai and 4> = 4>{i + 6(j) with the index denoting the background field solution. 
Note also that the operators D]^. and D]^ denote the differential operator one obtained from the linearized perturbation 
equation with all fields evaluated at the background solutions. The exact form of these differential operators will be 
shown later in the following arguments. 

One is looking for stability conditions that the field parameters must obey in order to keep the evolutionary 
solution stable. One can show that class (a) and class (b) solutions are good enough to cover all domain of stability 
conditions. Let us denote the domain of solutions to class (a), (b) and (a+b) stability conditions as 5(a), S'(6), 
and S{a + h) respectively. Specifically, the definition of these domains are defined by S{a) = {5ai\D]^.5aj = 0}, 
S(b) = {5(j)\D'^5(l) = 0}, and S{a + h) = {{5ai,5(l))\DlMj + D^^5(l) = 0}. 

Therefore, one only needs to show that S{a) n S{h) C S{a + h). This is because that "Dl^ Suj — and Dl^Scj) — 0" 
imply that ''Dl^Saj + D'^S<j) = 0". On the other hand, ""Dl^Saj + D'^S(j)h = 0" does not imply that ''D^Mj = 
or D]^54) — 0". Hence class (a) and class (b) solutions cover all the required stability conditions wc are looking for. 
Hence we only need to consider these two separate cases for simplicity. 

In addition, one knows that any small time-dependent perturbation against the metric field is known to be equivalent 
to a gauge choice |2^. This can be clarified as follows. Indeed, one can show that any small coordinate change of the 
form x'^ = — will induce a gauge transformation on the metric field according to g'^^^ = g^n + D^e^j + -Dj^e^. 
Therefore, a small metric perturbation against a background metric is amount to a gauge transformation of the form 
a[ = ai + e^cLi for the Bianchi type I metric with = (constant, ei(i)). This is then equivalent to small metric 
perturbations. If a background solution is stable against small perturbation with respect to small field perturbations, 
one in fact did nothing but a field redefinition. 

If the background solution is, however, unstable against small perturbations, e.g. the small perturbation will grow 
exponentially as we will show momentarily, the resulting large perturbations can not be classified as small gauge 
transformation any more. Therefore, the stability analysis performed in the literatures pO|-p^ for various models 
against unstable background solution served as a very simple method to check if the system support a stable metric 
field background. This is the reason why we still perform a perturbation on the metric field for stability analysis even 
small perturbation is equivalent to a gauge redefinition. 

Note that one should also consider a more general perturbation with space perturbation included. The formulation 
is, however, much more complicate than the one we will show in this paper. We will focus on the time-dependent 
case for simplicity in this paper. The space-dependent perturbation analysis is still under investigation. The time- 
dependent analysis alone will, however, brings us many useful information for the stability conditions about the model 
we are interested. For example, we will show in the following subsection that the solution found in Ref. |^ remains 
stable as long as the scalar field falls close to any local minimum of the potential U{4>). Note again that the solution 
found in Ref. ^ is an exact solution only when U = constant. 

A. Dilaton Field Perturbation 

In this subsection, we will consider the perturbation S(j) with respect to the exact background solution c/jb given in 
the previous section. 

We will consider a general symmetry-breaking potential C/(0) which has at least one local minimum. The solution 
obtained in section II remains a good approximation as long as the dilation field is close to the local minimum (</> = v) 
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of the potential U such that Uq := U{v) = A. Therefore, a perturbation on (j) field can also be considered as a test 
to see whether such solution tends to stabilize the FRW space or not. We will show momentarily that the result 
indicates that FRW space tends to be a stable final state of any Bianchi type I space once the dilaton field approaches 
the local minimum of the symmetry-breaking potential. 

Due to the fluctuation Scj), the deviation of the dilaton potential, up to the first order, is 

Ui<j)B + H) U{^b) + d^Ui^s) <50 = A + d^U{(l)B) 5dp. (23) 
By keeping the metric fields unperturbed, one can show that the perturbation equations are 



= 0, (24) 



9^[/(0b)'50 = O, (25) 

^dt{VB S^) + He H + ^dpicj^B) Scf> = 0. (26) 
Vb 

Here we use the the Eq. (^l|) in above equations. In addition, once the dilaton fields falls close enough to one of 
the local minimum of the symmetry-breaking potential, the dynamics of the dilaton perturbation is inevitably much 
smaller than the scale factor perturbations. This confirm our claims earlier that class (a) perturbation and class (b) 
perturbation are not of the same order of magnitude. 

Note also that the Eq. ( p5| ) indicates that one of the necessary stable conditions for the background solution is that 
(j)B should be a local minimum of the dilaton potential U{<j)). In fact, close to local minimum condition is the only 
situation we are interested in this paper. 

Combining the equations of 5<j), one finds that the dilaton perturbation (50 is determined by 



dt(yB 



e40B 5^) + LvBe^'<'-dlU{cl>B) Scb = 0. (27) 



Taking the large time limit, i — > cxd, one can show that the background variables Vb and (/>b approach 

Vb ^ ye"*, cj^B ^ V, (28) 



where a := y/3A/2 and v := ln(2(^Q)/2. In addition, the equation (27) reduces to 



6(j) + a Scj) + f] 6(1)^0, (29) 

with f3 — d'pJ{(t>)\^=y/2n. For the case d'pJ{(t))\^=y = 0, i.e. /? = 0, the perturbations of dilaton and axion fields 
behave as (50 oc exp(— (at) and 5h cx exp(— 2(a;i) indicating both fields approach zero exponentially in the large time 

limit. On the other hand, the dilaton perturbation can be shown to be 54) (x exp(7i) with 7 — (^—a ± \J o? — 4/3^ /2 

for the case /3 7^ 0. This implies that \8<\)\ vanishes rapidly near the local minimum v when d'p{(4))\^=v > 0. As well, 
dh cx exp[(7 — a)t] also converge to zero exponentially. 

In short, the necessary condition for the stability of the background solutions is that the background dilaton field 
must be close to the local minimum v of the dilaton potential U{(j)). Consequently, the perturbed fields will approach 
zero exponentially provided that d'^U{(l))\^=v is non negative. This partially refiects the fact that the background 
solution remains a good approximated solution to the system as long as the dilaton potential U {(p) ~ constant near 

4> = V. 

Note that we know the evolutionary solution of the system only when the dilaton potential is a constant. Nonethe- 
less, the argument shown in this subsection indicates that the system tends to bring the universe to the FRW space 
once the dilaton field falls close to local minimum of the symmetry-breaking potential as t ^ 00. This result provides 
a convincing evidence for the formation of the fiat FRW space evolved from a Bianchi type I anisotropic space-time. 



B. Metric Perturbation 



We will study the stability of the background solution with respect to perturbations of the metric in this subsection. 
Perturbations will be considered for all three expansion factors via 
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tti ^ aBi + = aBi{l + 5bi), (30) 

hereafter. We will focus on the variables Sbi instead of Sui from now on for convenience. Therefore, the perturbations 
of the following quantities can be shown to be 

e^^eB^+Sk, e^0B + ^J2^^^^ ^Y.^B^ + ^J2^B^Sb,, V ^ Vb + VbJ^ ^b,. (31) 

i i i i i 

One can show that the metric perturbations Shi^ to the hnear order in Shi^ obey the following equations 

^ <55, + 2 ^ OBrSk - Se-"^" J2 = 0, (32) 

V p 

6b, + -^6b, + ySb,eB^ = 0, (33) 
Vb 



= 0. (34) 



Equation ( p4| ) indicates that J2i = 0. Therefore, Eq. (^2|) gives further that 

^(56,-0. (35) 

i 

Substituting it back to the other equations, one can show that 

5b, + -^5b, = 0. (36) 
Note that the background variables Vb and Obi approach 



Vb ^ ye"*, 0Br ^ VA/6, (37) 

in the large time limit as <: ^ oo. Here a := \/3A/2. Note also that we keep only leading order in the above equations. 
Therefore, the solution for Sbi can be found to be 

Sb, = Qe-"*, (38) 
together with the following constraint on integration constants Ct 

J2c,^0. (39) 

i 

Moreover, the asymptotic form, in the large time limit, of background scale factors is a^i a.^o exp[Q;i/3]/4-'^/^. 
Therefore, the "actual" fluctuations for each expansion factor, Sai = UBiSbi is 

- ?f (40) 

Hence Sai approaches zero exponentially since a is definite positive. Consequently, the background solution is stable 
against the perturbation of the graviton field. 

Note that the perturbations of the three different expansion factors in Bianchi type I cosmology are not independent 
due to the constraint (|3^). The constraint ( ^9|) indicates that scale factor perturbations along arbitrary directions 
are confined by the field equations. This also signifies the symmetry of the coordinate transformation among the 
scale factors. In fact, there should be only two independent coordinate perturbations reflecting the symmetry of the 
coordinate choice. 
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IV. CONCLUSION 



In summary, we investigate the stability condition for an anisotropic inflationary string solution which evolves to 
the FRW space. This solution is consistent with the significant requirements constrained by present astronomical 
observations — homogeneity, isotropy and accelerating expansion. We analyze in details the perturbed equations 
with respect to the dilaton fields and separately with respect to the expanding scale factors of the Bianchi type I 
space. Our result indicates that the cosmological model considered in this paper is stable for both cases under certain 
constraints. 

In the former case, the necessary condition for the stability of the background solutions is that the background 
dilaton field must be close to the local minimum v of the dilaton potential U{(j)). Consequently, the perturbed fields 
will approach zero exponentially provided that d'^U{(j>)\ci,=v is non negative. This partially refiects the fact that the 
background solution remains a good approximation to the system only when the dilaton potential U {(j)) ~ constant 
near (p = v. 

Note that we do not known the exact evolutionary solution of the system when the dilaton potential is not close 
to a constant. Nonetheless, we show that the system tends to bring the universe to the FRW space once the dilaton 
field falls close to local minimum of the symmetry-breaking potential as t — > oo. This result provides a convincing 
evidence for the formation of the flat FRW space evolved from a Bianchi type I anisotropic space-time which will 
work for models with a reasoable symmetry-breaking potential. 

For the later case, the perturbed expanding scale factors are also shown to approach zero exponentially in the large 
time limit. We also show that the perturbations of the three expanding scale factors are constrained by the field 
equations. 
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